Exchange interactions in (ZnMn)Se 
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One of the remarkable properties of the II- VI diluted magnetic semiconductor (ZnMn)Se is the 
giant spin splitting of the valence band states under application of the magnetic field (giant Zeeman 
splitting). This splitting reveals strong exchange interaction between Mn moments and semicon- 
ductor states. On the other hand, no magnetic phase transition has been observed for systems 
with small Mn content up to very low temperatures. The latter property shows weakness of the 
exchange interaction between Mn moments. In this paper, the local density approximation (LDA) 
and the LDA-I-U techniques are employed to study exchange interactions in (ZnMn)Se. Supercell 
and frozen-magnon approaches applied earlier to III-V diluted magnetic semiconductors are used. 
It is found that both LDA and LDA-I-U describe successfully the combination of the strong Zeeman 
splitting and weak interatomic exchange. However, the physical pictures provided by two techniques 
differ strongly. A detailed analysis shows that the LDA-I-U method provides the description of the 
system which is much closer to the experimental data. 
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I. INTRODUCTION 

The perspective of using the spin of electrons in the 
semiconductor devices promises to revolutionize mod- 
ern electronicsi A necessary component of a spintronic 
process is the spin-injection into semiconductor. This 
demand created a need for ferromagnetic materials on 
the semiconductor basis with strong spin-polarization of 
the carriers and high Curie temperature. After recent 
discovery^ of ferromagnetism in (GaMn)As with Curie 
temperature as high as 110 K much attention is de- 
voted to the study of the diluted III-V semiconductors 
as possible sources of the spin-polarized electrons. Also 
the interest to more traditional II- VI semiconductors has 
been revived since, first, the theoretical studies shaw that 
the II- VI systems possess the potential for high Curie 
temperatureSii, second, the study of the II- VI diluted 
magnetic semiconductors (DMS) deepens understanding 
of the exchange interactions in other types of DMS^*^ 
and, third, the II- VI systems, in particular (ZnMn)Se, 
are used in spin-injection experiments as a source of spin- 
polarized charge carriers liiSiS 

In (ZnMn)Se with low Mn concentration no magnetic 
ordering has been experimentally detected up to very 
low temperaturesi^ However, the application of mag- 
netic field leads to the observation of the so-called giant 
Zeeman splitting for the states of the semiconductor ma- 
trix. Because of this large spin splitting (ZnMn)Se is an 
efficient source of highly polarized carriersj^i^ 

A commonly accepted explanation of the giant Zeeman 
splitting in (ZnMn)Se relies on the following picture. In 
the system there is a strong exchange interaction between 
the Mn 3d states and the states of the semiconductor. 
This interaction does not, however, lead to the ordering 
of the Mn moments. The application of a magnetic field 
aligns the Mn moments and results in the observation of 
a giant spin-splitting. The term 'giant' arises here from 
the comparison of two energy scales. The characteristic 
energy of the magnetostatic interaction of the magnetic 
field of 1 T with the spin of \^jlb amounts to 0.004 mRy 



and is up to five orders of magnitude smaller than the 
splittings detected experimentally.^' The coexistence of 
a very strong exchange between the Mn 3d and semicon- 
ductor states and a very low temperature of the magnetic 
phase transition makes the II- VI DMS an interesting lab- 
oratory for studying the physics of exchange interactions. 

Much efforts has been paid to the theoretical studies 
of the DMS of the II- VI type (see, e.g., reviewo'^i^^i^"). 
Most of these studies are based on a model-Hamiltonian 
approach. [See, e.g., Refiii for the model of bound mag- 
netic polarons and Refs&i^ for the Zcner model. The lat- 
ter can be considered as a continuous- medium limit of the 
well-known Rudermann-Kittel-Kasuya-Yosida (RKKY) 
approach. Spalek et al^'^ used Anderson's approach to 
superexchange and have shown that few adjustable pa- 
rameters of the theory are sufficient to describe the in- 
terplay between different types of exchange interactions 
in the system. Reach experimental information (e.g^ in 
the case of (ZnMn)Se) is helpful in the selection of the 
values of the parameters.] 

The developments in the methods of the density func- 
tional theory (DFT) accompanied by fast increasing com- 
puter power allow now for parameter-free calculation of 
the electronic properties of very complex systems. In the 
case of II- VI DMS the number of the DFT studies of the 
exchange interactions is still very restricted (for excep- 
tions see, e.g., RefsAi^ii^). Most of the calculations have 
been performed with the use of a virtual crystal approxi- 
mation or a single-site coherent-potential approximation 
(CPA)iJ4 These schemes are convenient and efficient in 
the investigation of the systems with varying concen- 
tration of impurities. They, however, do not take into 
account the atomic short-range order and the increas- 
ing distance between impurities with decreasing impu- 
rity concentration. Therefore it is important to combine 
the calculations within the virtual-crystal and CPA tech- 
niques with the studies making more detailed account for 
the positions of the impurity atoms. Such studies can be 
performed with the use of large supercells of the semi- 
conductor crystals liiiSiiSiSS The aim of this paper is the 
investigation of the exchange interactions in (ZnMn)Se 
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on the basis of the supercell approach. 

An important question concerns the role of the in- 
traatomic correlations in the Mn 3d states. Within the 
model-Hamiltonian approaches the 3d states are usually 
considered as atomic-like and strongly correlated. This 
treatment is rather different from the treatment within 
LDA. To study the role of the electron correlations in 
the Mn3d shell we use the LDA+U method — designed 
to take into account the on-site Coulomb interaction U. 

The purpose of this paper is two-fold. On the one hand, 
by detailed DFT calculation of the electronic properties 
of (ZnMn)Se we aim to provide deeper insight into the 
physics of the system. On the other hand, by compari- 
son with experimental data of the calculational results 
we aim to draw the conclusion which of the two ap- 
proaches, LDA or LDA-I-U, provides better description 
of (ZnMn)Se. Both components of the purpose are of 
strong importance for future studies of DMS within the 
DFT. 



FIG. 1: Two scenarios of the relation between the net mag- 
netization due to the localized moments and the exchange 
splitting of the valence state. In the mean-field (Stoner-like) 
scenario the splitting is proportional to the net magnetiza- 
tion. The distance between wave lines shows schematically 
the value of the exchange splitting. In the non-mean-field 
(non-Stoner) scenario the spin of the valence electrons fol- 
lows locally the direction of the Mn moments. In this case 
the exchange splitting is present also in the case of zero net 
magnetization. 
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II. CALCULATIONAL APPROACH 

In the calculations we follow the scheme described in 
Refii^. This scheme is based on the supersell approach 
where one of the Zn atoms in a supercell of zinc-blende 
ZnSe is replaced by the Mn atom. The calculations are 
performed for four values of the concentration x: 0.25, 
0.125, 0.0625, and 0.03125. 

The calculations were carried out with the augmented 
spherical wavesS^ (ASW) method within the LDA and 
LDA-I-U approaches. In all calculations the lattice pa- 
rameter was chosen to be equal to the experimental lat- 
tice parameter of ZnSe. Two empty spheres per formula 
unit have been used in the calculations. The positions of 
empty spheres are (0.5, 0.5, 0.5) and (0.75, 0.75, 0.75). 
Radii of all atomic spheres were chosen to be equal. De- 
pending on the concentration of Mn, the super cell is 
cubic (x=25%, ax ax a, and x=3.125%, 2a x 2a x 2a) or 
tetragonal (x=12.5%, ax ax2a and 6.25%, 2a x 2a x a). 

The LDA-I-U calculations were performed with U ~ 
0.3RyM 



A. Spin-projected densities of states 

The spin splitting of the valence-band states is a result 
of the interaction between these states and the Mn3d 
states. The directions of the Mn moments are disordered 
in the absence of magnetic field and become increasingly 
ordered with increasing value of the applied field. There 
are two possible scenarios for the relation between the 
ordering of the Mn moments and the spin splitting of the 
valence band states (Fig. First scenario is a mean- 
field (Stoner-like) type of the relation. In this case the 
valence-band states experience an average exchange field 
of the Mn moments that is proportional to the net mag- 
netization in the Mn subsystem. Complete disordering 
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of the Mn moments leads to the disappearance of the net 
magnetization and, respectively, of the spin splitting. In 
the second scenario, the spin of the valence states follows 
locally the spins of the Mn atoms. Therefore, the spin po- 
larization and exchange splitting do not disappear with 
disappearance of the net magnetization. An experimen- 
tal example of the non-Stoner behavior of the exchange 
splitting is discussed by KiskerfA for the case of iron with 
thermally disordered atomic moments. In (ZnMn)Se the 
experimental data is treated in favor of the mean-field 
scenario. 

To study both scenarios within the DFT we will need 
to calculate densities of states (DOS) projected on differ- 
ent spin-quantization axes. The calculation is performed 
as follows. The wave function of a given electron state is 

V'i(r) 



considered to be two-component spinor 



V'2(r) 



The 



spin components ipilr) and "02(1") are written with re- 
spect to a chosen axis. The spin-quantization axis does 
not change within atomic spheres but can vary from atom 
to atom. The integral 



(1) 



gives the part of the state Vik corresponding to atom 
1/ and spin-projection s. Here k is the wave vector, i 
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numbers the energy bands. The integration is carried 
out over the z/th atomic sphere. The wave functions are 
normahzed in the unit cell: X^si/ '^fkiy — 1- The partial 
DOS for given v and s is given by formula 
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BZ 



5{e - e,k) (2) 
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where ftsz is the volume of the Brillouin zone (BZ). To 
calculate the partial DOS with respect to another quan- 
tization axis the electron wave functions are subjected to 
the transformatior 



^i(r) 
-02 (r) 



r e 



(3) 



where Ui, is the spin-i transformation matrix corre- 
sponding to the rotation of the axis of the vth atom. 
Next, the procedure defined by Eqs. H1I2|) is performed 
for the spinor components of the transformed functions 



B. Interatomic exchange parameters and Curie 
temperature 

To describe the exchange interactions between Mn mo- 
ments we use an effective Heisenberg Hamiltonian of clas- 
sical spins 



Heff - - 



(4) 



where Jij is an exchange interaction between two Mn 
sites (i,j) and is the unit vector pointing in the direc- 
tion of the magnetic moment at site i. 

To estimate the parameters of the Mn-Mn exchange in- 
teraction we perform calculation for the following frozen- 
configurations: 



magnor 



0j = const, 



q 



(5) 



where 6i and 0i are the polar and azimuthal angles of 
vector Bi, Hi is the position of the ith Mn atom. The 
directions of the induced moments in the atomic spheres 
of Ga and As and in the empty spheres were kept to be 
parallel to the z axis. 

It can be shown that within the Heisenberg model 
the energy of such configurations can be represented in 
the form 



(6) 



where Eq does not depend on q and J(q) is the Fourier 
transform of the parameters of the exchange interaction 
between pairs of Mn atoms: 



>/(q) = X! exp(iq • Roj)- 



(7) 



In Eq. © angle 9 is assumed to be small. Using J(q) one 
can estimate the energies of the spin- wave excitations^^: 



where M is the atomic moment of the Mn atom. Per- 
forming back Fourier transformation we obtain the pa- 
rameters of the exchange interaction between Mn atoms: 



(9) 



The Curie temperature was estimated in the mean-field 
(MF) approximation 
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(10) 



We use rigid band approach to calculate the exchange 
parameters and Curie temperature for different electron 
occupations. We assume that the electron structure cal- 
culated for a DMS with a given concentration of the 3d 
impurity is basically preserved in the presence of defects. 
The main difference is in the occupation of the bands 
and, respectively, in the position of the Fermi level. 



III. CALCULATIONAL RESULTS 
A. Densities of states and exchange splittings 
1. LDA 

We begin the discussion of the calculational results 
with consideration of the DOS of the ferromagnetic 
(ZnMn)Se (Fig. EJ. Compared with pure ZnSe, the 
replacement of a Zn atom by a Mn atom adds five 3d 
spin-up energy bands to the valence band of the system. 
Since the number of Mn 3d electrons is also five no carri- 
ers appear either in the valence band or in the conduction 
band. 

The values of the calculated spin moments are collected 
in Table HI The moment in the Mn sphere and the in- 
duced moment on the neighboring Se atom are practically 
independent of the Mn concentration. The moment per 
supercell is exactly 5/is since there are extra five filled 
spin-up bands compared with the nonmagnetic ZnSe. 

At the top of the valence band there is strong negative 
exchange splitting: the upper edge of the spin-up DOS 
lies higher in energy than the upper edge of the spin- 
down DOS ( Fig. The spin splitting increases with 
increasing Mn concentration. The origin of this splitting 
can be understood from the analysis of the partial Mn3d 
DOS (Fig. ISJ. Indeed, a higher energy of the spin-up 
valence-band edge results from the contribution of the 
Mn3d spin-up states. On the other hand, there is no 
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TABLE I: Magnetic moments in Zni-ajMn^jSe. There are 
shown the Mn moment, the induced moment on the nearest 
Se atoms, and the magnetic moment of the super cell. All 
moments are in units of /is. 
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FIG. 2: The DOS of Zni_^Mn^Se. The DOS is given per 
unit cell of the zinc-blende crystal structure. The DOS 
above(below) the abscissas axis corresponds to the spin- 
up(down) states. 



(ZnMn)Se 




strong contribution of the Mn spin-down states to the 
valence band. Unoccupied spin-down Mn3d bands lie 
in the semiconducting gap, close to the bottom of the 
conduction band. 

The exchange splitting of the states at the top of the 
valence band of the ferromagnetic (ZnMn)Se is in qual- 
itative agreement with the observation of the giant spin 



FIG. 3: The partial Mn3d-D0S for Zm-^^^Mn^i^Se. The DOS 
is given per Mn atom. 
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splitting under the appUcation of the magnetic field. 

To understand the dependence of the electron struc- 
ture on the magnetic configuration we performed cal- 
culation for the antiferromagnetic configuration of the 
Mn moments for x = 3.125%. The antiferromagnetic 
structure is characterized by the largest angle between 
the neighboring Mn moments and, therefore, is the state 
most different from the ferromagnetic one. The com- 
parison of the DOS of the ferromagnetic and antiferro- 
magnetic configurations for x = 3.125% (Figs. [21 and 
lU shows that they are similar. The main difference in 
the antiferromagnetic DOS is the appearance of a small 
energy gap which separates the upper part of the valence- 
band states from the rest of the valence band. The spin- 
splitting obtained in the ferromagnetic DOS is present 
also in the antiferromagnetic DOS as the splitting be- 
tween the top of the separated (impurity) band and the 
top of the valence band. Figures 13 and ^ show that the 
partial DOS of the Mn atoms are only weakly dependent 
on the magnetic configuration and preserve fully the lo- 
cal spin splitting. Obviously the main features of the 
LDA-DOS cannot be treated in terms of the Stoner-like 
mean- field picture (Fig. 

To study in more details the relation between the elec- 
tron structure and magnetic configuration we perform 
the analysis of the partial DOS of the atoms of the semi- 
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FIG. 4: Total DOS and partial Mn3d-DOS for the antiferro- 
magnetic Zni_2,Mn2,Se with x — 3.125%. The spin-projected 
DOS are calculated with respect to the local atomic quanti- 
zation axes. 
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FIG. 5: Fragment of the partial DOS of the Se-I atom at the 
top of the valence band for Zno.g6875Mno.o3i25Se. Calculations 
are performed for seven magnetic configurations with different 
net magnetization. The upper part shows the spin-projection 
on the global z axis. The lower part gives the spin-projections 
on the axis parallel to the direction of the magnetic moment 
of the nearest Mn atom. No correlation between the spin- 
splitting and the net magnetization can be established. 
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FIG. 6: Fragment of the partial DOS of the Se-II atom. See 
the caption of Fig. |5] for details. Short vertical lines in the 
upper panel show the center of gravity for the correspond- 
ing spin-projection. The decrease of the spin-splitting with 
decreasing net magnetization can be established. 



conductor matrix. As we show below, the behavior of the 
DOS of different atoms ranges from a highly non-Stoner 
one to the behavior well described by the mean-field pic- 
ture. 

In Figs. I5I5I we present the partial DOS for two dif- 
ferent Se atoms. The DOS is calculated for magnetic 
configurations depicted schematically in Fig. |7| where 
three successive Mn atoms along the z axis are shown. 
The angles between neighboring Mn moments vary from 
to 180° with a step of 30°. Correspondingly, the net 
magnetization vary from maximal to zero. 

First of the two Se atoms, Se-I, is at position {j, |, |) 
and is the nearest neighbor of the Mn impurity situ- 
ated at (0,0,0). The second atom, Se-II, is at position 
(X' X' x) ^^'^ belongs to the 4th coordination sphere 
of Se atoms. (To remind, the supercell in this case is a 
cube of the size 2a.) The spin-projected DOS is presented 
with respect to two different quantization axes. The first 
(local) axis is parallel to the direction of the magnetic 
moment of the nearest Mn atom. The second (global) 
axis is directed along the net magnetization. The DOS 
projected on the local atomic quantization axis provides, 
in most cases, better insight into the physics of the sys- 
tem. However, in the experiments determining the spin- 
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splittings the spin-quantization axis is usually the global 
one. 

The DOS shown in Figs. I5I6I reveal strong difference 
between two Se atoms. The spin of the electron states 
of the Se-I atom follows almost perfectly the spin of the 
neighboring Mn atom (Fig. (S)). This is evidenced by 



6 



FIG. 7: Schematic picture of the calculated magnetic 
configurations. Angle 6 assumed the following values: 
0°, 15° ,30° ,45° ,60°, 75° ,90°. Calculations were performed for 
(ZnMn)Se with the Mn concentration of a; = 3.125%. 
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FIG. 8: The exchange splitting calculated from the partial 
DOS of the Se-II atom as a function of the net magnetization. 
Negative value of the splitting reflects higher energy position 
of the spin-up states. The dashed straight line is a guide for 
the eye. 
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the negligibly small local spin-down DOS. The smallness 
of the local spin-down DOS holds up to the largest angle 
between Mn moments. Considered from the viewpoint of 
the global quantization axis, both spin-up and spin-down 
DOS of Se-I have similar shape but different weights (am- 
plitudes). The relative weight of the spin-down DOS in- 
creases from zero for = 0° to 1 for 6* = 90°. Because of 
the similarity of the form of the global spin-up and spin- 
down DOS the variation of the DOS with the change 
of the net magnetization cannot be treated in terms of 
the varied spin-splitting. The changes in the global DOS 
take the form of the redistribution of the weight between 
the spin-up and spin-down DOS. This behavior is prin- 
cipally different from the behavior expected within the 
mean-field picture. 

On the other hand, for Se-II we get a strong depen- 
dence of the local DOS on the magnetic configuration of 
the Mn moments (Fig. With increasing angle between 
Mn moments the contribution of the spin-down DOS in- 
creases. This happens because the spin of the states of 
Se-II deviates increasingly from the spin of the nearest 
Mn atom responding to the influence of other Mn atoms. 
Now, the shape of the spin-up and spin-down DOS calcu- 
lated with respect to the global quantization axis differ 
strongly. To characterize this difference in terms of the 
spin-splitting we calculated the centers of gravity for both 
spin-DOS in the energy region at the top of the valence 
band (Fig. IS]). In a good approximation the spin-splitting 
is proportional to the net magnetization (Fig. IS)). There- 
fore the properties of the electron structure of Se-II can 
be described within the mean-field picture. 



FIG. 9: The exchange splitting at the top of the valence 
band (F point). Calculations are performed within LDA and 
LDA-I-U. The solid line corresponds to the experimental value 
of the Jpd parameter.— 
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Since the relation between the exchange splitting and 
net magnetization varies from atom to atom, in exper- 
iments probing different parts of the electron structure 
this relation can appear different. 

Important role in the experimental determination of 
the strength of the exchange interaction between the Mn 
moments and valence band states in DMS is played by 
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the magneto-optical measurements of the spin splitting 
at the top of the valence band (F point of the Brillouin 
zone)»ifi In Fig. IHlwe show the calculated dependence of 
the splitting on Mn concentration. The LDA results for 
the exchange splitting at the F point disagree with exper- 
iment in two respects. First, the mean-field picture which 
assumes proportionality between the exchange splitting 
and the net magnetization in the system does not apply 
to the LDA results. The LDA-splitting is well described 
by a linear function with a finite value in the limit of 
X — > 0. Such form of the dependence is the result of 
the presence of the Mn3d states at the Fermi level. Sec- 
ond, the exchange splitting is substantially larger than 
the splitting obtained experimentally. Since the LDA re- 
sults cannot be described within the mean-field picture 
the use of the formula 

J,. ^ (11) 

gives values of the Jpd parameter that depend on con- 
centration X. For illustration, the value of 3^^ obtained 
for X = 0.125 according to Eq. is about two times 
larger than the experimental value. In Eq. Hll|l . /S.E is 
the exchange splitting and S is the atomic spin moment 
of Mn. 

2. LDA + U 

Introducing of the Hubbard-U into calculational 
scheme results in a strong shift of the Mn3d spin-up 
states to lower energies (Fig. I1UI11|I . At the top of 
the valence band there is still an admixture of the Mn3d 
states. It is, however, very weak, especially for the spin- 
up states. Nevertheless, one can notice the hybridiza- 
tional repulsion of the valence-band spin-down states 
from the spin-down Mn3d states lying in the semicon- 
ductor gap. This repulsion is an important factor lead- 
ing, in agreement with experiment, to negative exchange 
splitting at the top of the valence band. Note that in 
the LDA calculation the negative exchange splitting has 
different origin: the presence of the spin-up Mn3d states 
at the top of the valence band. 

To study the relation between net magnetization and 
exchange splitting we performed calculation for the an- 
tiferromagnetic configuration of the Mn moments for 
X = 3.125% (Fig. 1121). In contrast to LDA (Fig. gj), 
no exchange splitting of the states at the top of the 
valence band is obtained in the antiferromagnetic case. 
This property is in agreement with the mean-field pic- 
ture (Fig. To investigate this property further we 
performed the LDA+U calculation of the exchange split- 
ting at the F point of the BZ as a function of the Mn 
concentration x (Fig. O . The result obtained within the 
LDA-I-U scheme is in very good agreement with exper- 
iment concerning both the mean-field character of the 
dependence and the magnitude of the exchange parame- 
ter Jpd- Summarizing the study of the exchange interac- 
tion between Mn3d and valence band states we draw the 



FIG. 10: The DOS of Ztii-^^Mn^^Se calculated within the 
LDA+U approach with U=0.3Ry. (To compare with Fig. |21 
presenting LDA-calculation.) 
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conclusion that the LDA+U approach provides for this 
property much better agreement with experiment than 
the LDA approach. 



B. Interatomic exchange interactions 

Now we turn to the discussion of the exchange inter- 
actions between Mn moments and address the question 
why large Mn moments and strong p-d exchange do not 
result in the case of (ZnMn)Se with low Mn content in 
a sizable magnetic phase-transition temperature. To get 
deeper insight into formation of the interatomic exchange 
interactions in the system we performed calculations for 
different band occupations (Fig. I13|l . The number of 
electrons varied from n — —2 (two electrons per super- 
cell less) to n = 0. Negative values of the Curie tem- 
perature in Fig. 1131 indicate an instability of the ferro- 
magnetic state due to dominating antiferromagnetic in- 
teractions. The calculations have been performed within 
both LDA and LDA+U approaches. In Fig. E]we show, 
for X — 3.125%, the dependence of the main interatomic 
exchange parameters as a function of the valence-band 
occupation. 
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FIG. 11: The partial Mn3d-D0S for Zni_^Mn^Se calculated 
within the LDA+U approach with U=0.3Ry. (To compare 
with Fig. 13 presenting LDA-calculation. Notice the change 
in the scale of the ordinate axis that reflects strong change in 
the DOS.) 



FIG. 12: Total DOS and partial Mn3d-D0S for the antifer- 
romagnetic Zni_a;Mna;Se with x — 3.125% calculated within 
LDA+U approach. The spin-projected DOS are calculated 
with respect to the local atomic quantization axes. 
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We begin with the discussion of the features common 
for both LDA and LDA+U calculations. Analysis of the 
calculated ( Fig. I13|l shows that in the case of com- 

pletely filled valence band and empty conduction band 
(n = 0) the main exchange interactions are antiferromag- 
netic. The antiferromagnetic character of the interaction 
agrees with the commonly excepted picture that the in- 
teratomic exchange interaction between magnetic atoms 
in an insulating system is dominated by the antiferro- 
magnetic superexchange. This interaction is considered 
to be mediated by the states of the intermediate non- 
magnetic atoma^ or the states of the completely filled 
bands 

Remarkable, however, is very small value of the ex- 
change interactions for n — 0. With decreasing number 
of the holes the absolute value of all interatomic exchange 
interactions becomes very small. This result is in good 
correlation with the failure to experimentally detect the 
spin-ordering up to very low temperatures^ It is also in 
agreement with the perturbative calculation by Larsen et 
al for Cdi_:^Mn^Tei^ 

The calculations show that the appearance of holes 
results in increasing ferromagnetic interactions. This is 



FIG. 13: T^^^ of (ZnMn)Se as a function of the electron 
number n. n = corresponds to the nominal electron number 
in (ZnMn)Se. For n = there is no charge carriers in the 
system. 
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FIG. 14: The main interatomic exchange parameters of 
(ZnMn)Se with Mn concentration of a; = 3.125% as a function 
of the electron number. 
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reflected by the property that the minimum of the esti- 
mated Tc is at n = (Fig. I13|l and correlates with ex- 
perimental observation of the ferromagnetism with small 
Curie temperature in p-doped (ZnMn)Te.^^ 

Although both LDA and LDA-I-U give very weak neg- 
ative exchange interactions for completely filled valence 
band the form of the dependence of the interatomic ex- 
change parameters and correspondingly Curie temper- 
ature on the number of holes differs strongly for these 
two techniques (Figs. I13I14|I . Comparison of the cal- 
culated exchange parameters shows that the main dif- 
ference between LDA and LDA-I-U concerns parameter 
Jo22 that describes the exchange interaction between the 
Mn atoms separated by vector (022): the LDA predicts 
much stronger increase of the ferromagnetic interaction 
than LDA-I-U. The strength of the ferromagnetic inter- 
actions obtained within LDA seems to be strongly over- 
estimated since the large values of the LDA-Tc (Fig. I13() 
does not correlate with the experimental data available. 
On this basis, we draw the conclusion that the LDA-fU 
scheme gives better description of the character of the 
dependence of the interatomic exchange interactions on 
the number of holes. 

This last conclusion might seem to be expected since 
in Sec. IIII Al we have seen that the value of the Jpd pa- 
rameter playing an important role in the mediation of the 
ferromagnetism is overestimated by LDA. The situation 



is, however, more complex since a larger Jpd parameter is 
accompanied, in the case of LDA, by a stronger spatial lo- 
calization of the valence-band hole states about the Mn 
atom. Increased localization of the holes produces the 
trend to decreasing Tc- Which of two factors, Jpd or 
hole-localization, prevails depends on the peculiar inter- 
play of the details of the electron structure of the specific 
system studied. The calculations for various III-V DMS 
show that the account for Hubbard-U can lead to both 
increase and decrease of the Curie temperaturen^Sr 

To complete the consideration of the exchange inter- 
actions in (ZnMn)Se we discuss the physical reason for 
a very weak interatomic exchange interactions for un- 
doped (ZnMn)Se (n = in Fig. [TSJ. To remind, this 
result is common for both LDA and LDA-I-U and is very 
important for understanding of the failure to experimen- 
tally detect any magnetic phase transition up to very 
low temperatures. Note, that the behavior of the II- 
VI DMS is different from that obtained for the III-V 
DMS (GaMn)As. In (GaMn)As, the same calculational 
scheme gives-'^^ for the completely filled valence band and 
X = 3.125% interatomic antiferromagnetic interactions 
that are about 20 times larger that in (ZnMn)Sc (Fig. 

To understand the weakness of the calculated su- 
perexchange in (ZnMn)Se and the difference between 
(ZnMn)Sc and (GaMn)As we invoke a tight-binding 
model of noncoUinear magnetic configurations. We con- 
sider helical configurations of the atomic moments and 
study the dependence of the band energy of the system 
on magnetic structure. The helical structures are defined 
by the formula 

e„ — (cos(q • R„) sin 9, sin(q • R„) sin 9, cos 9) (12) 

where R„ are the lattice vectors, q is the wave vector of 
the helix, e„ is the unit vectors in the direction of the 
magnetic moment at site R„, polar angle 9 gives the de- 
viation of the moments from the z axis. The helical struc- 
tures allow to describe broad range of magnetic configu- 
rations from coUinear ferromagnetism (0 = or q = 0) 
to collinear antiferromagnetism (q = and 9 — 90°, 
K is a reciprocal lattice vector). 

The tight binding method for spiral structures was dis- 
cussed in its general form in Refi^. By neglecting the 
difference in the spatial dependence of the basis func- 
tions with opposite spin projections and by preserving 
only the single-center matrix elements of the exchange 
potential we arrive at the following simple form of the 
secular matrix i?(k)^: 



iJ(k)= 



cos2 f + sin^ § H+ 
-\fim9{H^ -H+ 



iA -\sm9{H_-H+) 

sin^ f if- + cos^ I if+ + iA 



(13) 



10 



where — Ho(k— ^q), iJ+ = Ho(\s.+ iq), and matrix 
Ho (k) describes spin-degenerate bands of a non-magnetic 
crystal; A is the diagonal matrix of the on-site exchange 
splittings. In LDA-I-U scheme A includes also ^. Secu- 
lar matrix describes a many-band system and takes 
into account the hybridization between the states of dif- 
ferent atoms. Through the hybridization between the 
Mn states and the states of the semiconductor matrix 
the spin polarization is transmitted to the nonmagnetic 
atoms. Note, that in contrast to the two-band model used 
in Ref.^^ the matrix (|13l) includes all relevant bands. In 
particular the Mn3d states are assumed to be included. 
This makes model (|13|l conceptually similar to our ASW 
calculations. 

The property of the secular matrix H13|) that is impor- 
tant for us reads 

/ dkSp[H{k)]=2 [ dkFo(k) (14) 
Jbz Jbz 

that is the trace of the matrix -ff (k) does not depend on 
the magnetic configuration. The variation of the mag- 
netic structure changes the energy of individual electron 
states. The changes of different states, however, compen- 
sate. In the case the trace of the secular matrix represents 
the total energy of the system the invariance with respect 
to the magnetic configuration means that all effective in- 
teratomic exchange interactions are negligible. 

The property given by Eq. (|14() applies to the total 
energy of a system if all bands described by the tight- 
binding secular matrix are occupied. This condition 
can be fulfilled only in the case that the hybridization 
between the occupied and empty states is weak and can 
be neglected. The weakness of the hybridization allows to 
include into the secular matrix the occupied states only. 

An attempt to use Eq. (|14|l for the explanation of the 
weakness of the superexchange in (ZnMn)Se leads imme- 
diately to the following difficulty. According to Eq. H13() 
both spin-up and spin-down Mn 3d states must be in- 
cluded into the secular matrix to fulfill Eq. 1)14(1 . Since 
only the spin-up Mn 3d states are occupied, the inclu- 
sion of the spin-down Mn 3d states, apparently, does not 
allow to relate the trace of the matrix to the energy of 
the system. In the case of (ZnMn)Se this difficulty can, 
however, be overcome if we notice that the spin-down 
Mn 3d states form very narrow energy bands lying in the 
semiconducting gap of ZnSe (Figs. I2I10|I . For example, 
the estimation for the LDA case shows that the states of 
these bands are strongly localized about the Mn atoms: 
more than 76% is located in the Mn spheres and more 
than 88% within the first coordination sphere of the Zn 
atoms. These states can be treated as evanescent states 
that are unable to mediate efficiently the exchange in- 
teraction between Mn atoms (Fig. I13|l . The evanescent 
character of the empty Mn 3d states allows to approxi- 
mately consider their contribution into the trace of the 
secular matrix [Eq. (|14l) ] as being independent of the 
magnetic configuration. This has as a consequence that 
the contribution of the occupied states is also approxi- 



mately independent of the magnetic configuration result- 
ing in weak effective interatomic exchange interactions. 

This consideration allows us to also explain the differ- 
ence between DMS on the GaAs and ZnSe bases. The 
semiconducting energy gap is substantially larger in the 
case of ZnSe. In the case of GaAs the influence of empty 
states is stronger and they must be included into the 
secular matrix H13|l . Therefore, the property described 
by Eq. H14|) does not apply to the occupied states that 
results in stronger superexchange. 

Some further comments are worth making here. First, 
the relation between the value of the semiconducting 
gap and the spatial extent of the exchange interactions 
has been many times discussed in the scientific literature 
within the framework of the perturbative treatment in- 
volving virtual transitions between the occupied states 
of the valence band and empty states of the conduction 
band (see. e.g., Refi^ for an early publication on this 
topic. See also an interesting comment on the relation 
between the range of the exchange interaction and an 
imaginary Fermi vector in recent paper by Pajda et al^). 
The model considered above captures basically the same 
physics by taking into account, in a non-perturbative 
manner, the hybridization between occupied and empty 
states. Only in the case the contribution of this hy- 
bridization into the response of the valence band states 
on the change of magnetic configuration is small the in- 
teratomic exchange interactions are weak. An important 
feature of model (|13|l is that it reflects the properties of 
the non-perturbative technique we used in the calcula- 
tion of the exchange interactions and, therefore, provides 
additional understanding of the calculational results. 



IV. CONCLUSIONS 

One of the remarkable properties of the II- VI diluted 
magnetic semiconductor (ZnMn)Se is the giant spin split- 
ting of the valence band states under application of the 
magnetic field (giant Zeeman splitting) . This splitting re- 
veals strong exchange interaction between Mn moments 
and semiconductor states. On the other hand, no mag- 
netic phase transition has been observed for systems with 
small Mn content up to very low temperatures. The lat- 
ter property shows weakness of the exchange interaction 
between Mn moments. In this paper, the local density 
approximation (LDA) and the LDA-I-U techniques are 
employed to study exchange interactions in (ZnMn)Se. 
Supercell and frozen-magnon approaches applied earlier 
to III-V diluted magnetic semiconductors are used. It 
is foimd that both LDA and LDA-I-U describe success- 
fully the combination of the strong Zeeman splitting and 
weak interatomic exchange. However, the physical pic- 
tures provided by two techniques differ strongly. A de- 
tailed analysis shows that LDA+U method provides the 
description of the system which is much closer to the 
experimental data. 
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